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ABSTRACT 

We present deep, large area B and r' imaging for a sample of 49 brightest cluster 
galaxies (BCGs). The clusters were selected by their x-ray luminosity and redshift to 
form two volume limited samples, one with (z) ~ 0.07 and one at (z) ~ 0.17. For each 
cluster the data cover 41' x 41'. We discuss our data reduction techniques in detail, 
and show that we can reliably measure the surface brightness at the levels of \xb ^29 
and /v ~ 28. For each galaxy we present the B and r' images together with the 
surface brightness profile, B — r' colour, eccentricity and position angle as a function 
of radius. 

We investigate the distribution of positional offsets between the optical centroid 
of the BCG and the centre of the X-ray emission, and conclude that the mass profiles 
are cuspy, and do not have extended cores. We also introduce a method to objectively 
identify the transition from BCG to extended envelope of intra-cluster light, using the 
Petrosian index as a function of radius. 

Key words: galaxies: fundamental parameters; galaxies: haloes; galaxies: photometry 



1 INTRODUCTION 

Brightest Cluster Galaxies (BCGs), form a unique popula- 
tion of objects both in their own right, and for the study of 
cluster formation and evolution. They are the most massive 
galaxies in the Universe, with typical masses of ~ 1O 13 M , 
comparable to that of a galaxy group, and luminosities of 
~ 10L* . BCGs are usually found close to the peak of cluster 
X-ray emission (e.g. Jones & Forman 1984) and at the centre 
of local density peaks (Beers & Geller 1983) , suggesting that 
BCGs are located at the bottom of their host cluster's po- 
tential well. BCG luminosity is also found to correlate with 
many global cluster properties such as X-ray temperature 
(e.g. Edge & Stewart 1991), and these facts taken together 
indicate that the origin of the BCG is closely linked to that 
of the cluster. 

The commonly favoured scenario for the formation of 
BCGs is through mergers of compact galaxy groups early in 
the history of the Universe (Merritt 1985). In this picture, 
BCGs form relatively quickly as the velocity dispersion in 
small groups is low allowing rapid merging. The simulations 
of Dubinski (1998) indicate that if this merging group falls 
into a collapsing cluster, the growing BCG will continue to 
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be fed by mass flows along filamentary large-scale structure 
as expected in hierarchical cosmological models. Some sup- 
port for this idea comes from observations that the major 
axes of BCGs are well aligned with both the X-ray isophotes 
of their host clusters and the cluster galaxy population (Bin- 
gelli 1982; Porter, Schneider & Hoessel 1991 and references 
therein). Zabludoff et al. (1993) report that many BCGs do 
not lie at the kinematic centres of their clusters, and con- 
clude that this is due to the BCG remaining at the kinematic 
centre of the sub-clump that formed it before or during its 
fall into the cluster. 

Many BCGs are found to have an extended, diffuse en- 
velope around them: these are classified as cD galaxies. cDs 
are usually found in the centres of aggregations of galaxies, 
strongly suggesting that the formation of the envelope is inti- 
mately tied to the group or cluster environment. Schombert 
(1988) argued that the excess light forms a component dis- 
tinct from the central galaxy and various authors have ar- 
gued that it arises from stars liberated by the stripping 
and/or disruption of cluster galaxies (Miller 1983; Malamuth 
& Richstone 1984). Linked with the problem of the nature of 
the cD envelope is that of intracluster light. This is thought 
to be due to a population of unbound stars pervading the 
cluster environment and orbiting in the cluster potential, 
yet its origin remains unclear. Often, the cD envelope and 
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ICL are not even separated in studies and we hope to in- 
vestigate whether this is a valid approach, or if meaningful 
distinctions can be made between the two components (see 
§5). 

Although the presence of ICL was detected as long ago 
as 1951 by Zwicky, it is only within the past decade or so 
that precise measurements have been made of its properties; 
this is partly due to its extremely low surface brightness, 
typically 1% of the night sky. It is observed to be a common 
component of clusters, contributing between 10% and 50% 
of the total optical cluster luminosity, and related to the 
dynamical state of the cluster (Arnabodi 2004; Willman et 
al. 2004; Feldmeier et al. 2004; Gonzalez et al. 2005; Zibetti 
et al. 2005). 

As noted by Zibetti et al. (2005), a large sample that 
allows generalisation of the properties of ICL is still lacking. 
Their work examines the ICL population at z = 0.2 — 0.3, 
which will be complemented by our lower redshift sample. 
The work of Gonzalez et al. (2005), consisting of a sample 
of 24 BCGs, has gone a major way to addressing this 
and Feldmeier et al. (2004b) are undertaking a survey of 
ICL in clusters not necessarily containing a cD galaxy. 
However, both studies use single pass-bands. One of the 
essential properties of ICL that needs to be pinned down is 
its colour, in order to shed more light on its origins. The 
simulations of Sommer-Larsen, Romeo & Portinari (2005) 
predict ICL colour to be similar to that of large ellipticals, 
whereas those of Willman et al. (2004) predict them to be 
bluer, though redder than dwarf galaxy populations. The 
current and forthcoming ICL studies of Krick, Bernstein 
& Pimbblet (2006), who have a sample of 10 clusters 
in V and r will help to address these issues, and we 
also hope to make headway in this area with our deep 
(fiB ~ 29 mag arcsec -2 ) imaging in the B and Sloan r' 
bands of a sample of 49 BCGs from a well-defined, X-ray 
selected sample of clusters. We intend this paper as an 
announcement of data release, though we are presenting 
images and radial profiles (these are available online at 
http:/ /www. nottingham.ac.uk/astronomy/research/bcg/ 
bcgdata.html) of our BCGs, and some basic analysis. In 
a further series of papers we will explore the properties of 
BCGs, cD envelopes, ICL and the relationships between 
these cluster components. 

Throughout this paper we assume a cosmology with 
Ho = 70km s _1 Mpc" 1 , D. = 0.27, and A = 0.73. 



2 DATA 

2.1 Sample Selection 

The main problem that has bedevilled work in the study of 
BCGs is choosing an appropriate sample. Using an optically 
selected sample, Aragon-Salamanca, Baugh & Kauffmann 
(1998) found a high rate of BCG evolution since z ~ 1, 
as expected in theoretical models, whereas the X-ray se- 
lected sample of Collins & Mann (1998) displays far less 
evolution. Burke, Collins & Mann (2000) showed that the 
two samples are consistent if there is a luminosity depen- 
dence on the rate of evolution, with only the BCGs of the 
lower Lx clusters showing significant mass evolution since 
z ~ 1. We have therefore selected our clusters from the well- 
defined Brightest Cluster Survey (BCS) of Ebeling et al. 
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Figure 1. The Lx — z distribution of the BCS visible during 
our observing run. The solid lines show the boundary of the high- 
luminosity volume-limited sample (0.05 < z < 0.2 and Lx > 
4.46 X 10 44 ergs — 1 ), and the dashed lines show the boundary of 
the low-redshift volume-limited sample (0.05 < z < 0.1 and Lx > 
1.1 X 10 44 ergs — 1 ). Large circles show clusters that we observed, 
with filled circles being those that we are presenting data for and 
open circles being those for which we could not process the images 
(see §2.2). Small points are clusters within our redshift and flux 
limits that we were not able to observe, and crosses show other 
clusters in the parent sample. 



(1998, 2000) as this will allow us to correct for luminosity- 
dependent effects. The BCS is an X-ray selected and X-ray 
flux limited set of clusters from the ROSAT All-Sky Survey, 
extending out to a z — 0.42. Our full sample of 60 clusters 
consists of a volume-limited subsample with 0.05 < z < 0.1 
and Lx > 1.1 x 10 44 ergs _1 , and a second volume-limited 
subsample with 0.05 < z < 0.2 and L x > 4.46 x 10 44 ergs _1 . 
The redshift range was specifically chosen to allow us to see 
the characteristic 'break' in the cD profile where the en- 
velope begins, and still allow us to perform accurate sky 
subtraction. This break is defined as the radius at which 
an excess of light is seen over the r 1//4 profile, and occurs 
at 50 - 500/i _1 kpc from the BCG centre (Schombert 1988). 
The Lx ~ z distribution of the BCS, along with our selected 
clusters (shown as filled and open circles), is shown in fig- 
ure Q Table 1 gives the X-ray properties of the BCG host 
clusters. 

X-ray selection guarantees that the clusters we are ob- 
serving are genuine massive bound systems, and not chance 
alignments of galaxies. This would affect examinations of 
BCG properties with, for example, cluster richness or the 
distribution of other cluster members. 

Another reason for using X-ray selected clusters is that 
data is available to allow us to ensure our BCG candidates 
are the central cluster galaxies we want to study. For a few 
clusters in our sample the BCG is not immediately obvious 
as the second-ranked cluster member appears to have a sim- 
ilar magnitude (Bautz- Morgan type II clusters). Knowledge 
of the position of the X-ray centroid therefore allows the 
most central cluster galaxy to be identified in these ambigu- 
ous cases. 
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Table 1. BCG/Cluster parameters. Redshifts and X-ray characteristics are taken from Ebeling et al. (1998, 2000). 



Cluster 


OtXRAY 1 


&XRAY 2 


otBCG 3 


Sbcg 4 


Offset (") 5 


Offset (kpc) 6 


L x (10 44 erg/s) 


z 


A0971 


10:19:55.2 


40:59:45.6 


10:19:52.1 


40:59:18.8 


54.2 


92.3 


1.44 


0.093 


A1126 


10:53:48.7 


16:50:31.2 


10:53:50.3 


16:51:03.0 


39.5 


62.6 


1.15 


0.086 


A1190 


11:11:28.6 


40:49:48.0 


11:11:43.6 


40:49:14.7 


228.2 


338.1 


1.47 


0.079 


A1246 


11:23:54.5 


21:29:16.8 


1 1 -23-58 8 


21:28:46.7 


71.3 


224.2 


7.62 


0.190 


A1589 


12:41:18.0 


18:33:03.6 


12:41:17.5 


18:34:28.2 


84.9 


114.8 


2.39 


0.072 
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57.7 
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4.96 


0.200 
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1 3 03-46 6 
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31.7 


38.2 


1.61 
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33-36 03 6 
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4.90 


0.188 


A1677 
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30-53 09 6 


1 3-05-50 8 


30-54-1 7 7 


75.8 


231.2 


6.24 


0.183 


A1728 


13:23:30.2 


11:17:45.6 


13:23:31.7 


11:18:08.1 


31.5 


52.9 


1.29 


0.091 


A1767 


1 3-36-07 7 


59 12-39 6 


1 3-36-08 2 


59:12:24.2 


17.5 


23.1 


2.47 


0.070 


A1775 


13:41:50.4 


26:22:55.2 


13:41:49.2 


26:22:24.1 


36.0 


49.0 


2.91 


0.072 


A1795 


1 3-48-52 3 


26-35-52 8 

J— yj .-)'). » i— . yj 


1 3-48-52 5 


26-35-35 3 


17.8 


21.1 


11.27 
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28:06:21.6 
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60.5 


84.9 
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29.6 


43.7 
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27-58-40 8 
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39.3 


45.8 
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0.061 
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6.2 
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2.40 
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18.39 
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A1927 
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25-37-40 8 
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52.5 


87.8 


2.14 


0.091 
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18:39:00.0 


14:54:31.5 


18:38:32.6 


28.6 


32.0 


1.38 


0.059 
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15-10-54 9 


5:43:12.0 


15-10-56 1 


5:44:41.8 


91.5 


131.2 


15.29 
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2.57 
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30-38-24 


15:21:20.6 
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3.95 
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27:42:51.9 


45.3 


61.7 


4.94 
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A2108 
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15-39-46 4 


1 7-50-09 4 
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1.97 


0.092 
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16:03:02.2 


25:24:14.4 
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0.089 
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29-53-31 2 


16:20:31.1 


29:53:28.2 
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0.097 
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69.0 


9.34 
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53.7 


80.3 


3.95 


0.080 
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17:17:45.9 
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78.2 
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A2457 


22:35:40.3 
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10:55:01.2 
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21:27:52.6 
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13:26:18.0 


0:13:33.6 


13:26:17.6 


0:13:17.9 


16.7 


25.6 


1.69 


0.082 


RXJ1442 


14:42:17.5 


22:18:03.6 


14:42:19.4 


22:18:11.5 


29.0 


51.4 


2.66 


0.097 


RXJ1750 


17:50:16.1 


35:04:58.8 


17:50:16.7 


35:04:57.9 


9.6 


27.8 


5.49 


0.171 


Z4905 


12:10:17.0 


5:23:31.2 


12:10:16.8 


5:23:09.6 


21.8 


31.4 


1.20 


0.077 


Z5029 


12:17:41.3 


3:39:32.4 


12:17:41.1 


3:39:21.7 


11.1 


15.6 


5.28 


0.075 


Z6718 


14:21:36.2 


49:32:38.4 


14:21:35.8 


49:33:03.0 


25.5 


34.1 


1.24 


0.071 


Z9077 


23:50:34.5 


29:31:51.6 


23:50:37.5 


29:29:07.6 


169.7 


295.7 


2.11 


0.095 



J R.A. of cluster X-ray emission peak 
2 Dec of cluster X-ray emission peak 
3 R.A. of BCG centre in r' band 
4 Dec of BCG centre in r' band 

Offset between X-ray peak and BCG centre in arcsec 
6 Offset between X-ray peak and BCG centre in kpc 
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2.2 Observations 

Our data were obtained during a single seven-night run in 
May/June 2003, using the 2.5m Isaac Newton Telescope 
(INT) at the Roque de Los Muchachos Observatory, La 
Palma, Spain. The images were taken with the Wide Field 
Camera (WFC), which is a four-CCD mosaic camera with 
a field of view of ~ 34' x 34', and a scale of 0".33/pixel. 
Each processed cluster image is composed of two offset WFC 
pointings; the BCG is centered on chip 4 of the camera in 
one pointing, and chip 2 in the other (see figure[5J . Chip 2 of 
pointing 1 does not form a contiguous part of the combined 
image and since each chip has a different gain, determining 
the correct sky level for chip 2 and still allowing for any 
large-scale sky gradients is extremely difficult. Hence, it is 
not included in the final mosaic. Rejection of this single chip 
gives a total field size of ~ 41' x 41'. Exposures were taken 
in Harris B (which closely mimics Johnson B — Salzer et al. 
2000) and Sloan r' filters, typically with all B mosaics be- 
ing a combination of 2 x 600s exposures, and the r' mosaics 
being 2 x 300s. Higher z clusters were observed for longer 
(see table 2) . Photometry was calibrated using Landolt stan- 
dard fields (Landolt 1992). The Cousins R magnitudes used 
by Landolt were converted to Sloan r' using the relations 
determined by Smith et al. (2002): 

' _ / V-0.44(B - V) +0.12 (V-R^1.00) , , 

r ~ \ V - 0.81(V - ii) + 0.13 (V-R>1.00) ( ' 

Conditions were photometric during the entire run, and 
seeing averaged at l".2. Sky surface brightnesses, large- 
scale surface brightness uncertainties and exposure times 
are given in table 2. The uncertainties give, in effect, the 
level to which our measured surface brightnesses are reli- 
able. A more complete explanation of the importance and 
calculation of these large-scale uncertainties is given in §3.4. 

Figure^shows that we observed all clusters in our lower 
z bin but that we are only ~ 60% complete in the higher 
z range. However, there is no obvious bias in the higher z 
selection. Out of the 60 clusters we did observe, we were not 
able to extract reliable surface brightness profiles for the 
BCGs for 11 of them (A1773, A2142, A2187, A2204, A2218, 
A2259, A2396, A2409, RXJ1720, RXJ1844 and Z8276) due 
to contamination by foreground sources. 



3 REDUCTIONS 

There are two critical factors which limit the measurement 
of low surface brightness features at large radii: the accuracy 
of flat-fielding and that of sky subtraction. The importance 
of these is pointed out in several other studies looking for 
intracluster light; for instance, Feldmeier et al. (2002) spent 
half their observing time taking blank sky images in order 
to construct sky flats, whilst Krick, Bernstein & Pimbblet 
(2006) used a third of their observing time for the same 
process. Due to the large number of observations we have, 
and the offsetting between pointings of the same cluster, we 
were able to construct our sky flats from our object images. 
This procedure is described in detail below, along with our 
sky subtraction algorithm. 



Table 2. Observational parameters. The left and right columns 
give the quantities for the B band and r' band observations, re- 
spectively. See §3.4 for more about large-scale errors (LSE). 



Cluster Sky level 1 LSE Exposure time 2 

B r' B r' B r' 



*Sky surface brightness and LSE in mags/arcsec 2 
2 Exposure time for a single pointing, in seconds 
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Figure 2. Mosaic combination of WFC pointings. Each of the 
final mosaics was made up of two pointings, one centred on chip 
4 of the WFC, the other centred on chip 2. The shaded (lower) 
pointing is pointing 2, the unshaded is pointing 1. Chip 2 of point- 
ing 1 is not included in the final mosaic as it does not overlap with 
any other part of the image. The numbers show the position of 
each of the WFC chips, which have active regions of 2048 X 4100 
pixels. The field of view of the mosaic is ~ 41' X 41'. 



3.1 Standard Reductions 

Data were reduced using a combination of IRAF 1 and Star- 
link 2 routines. Full two-dimensional debiassing was neces- 
sary, after which bad pixels were found (by median com- 
bining the frames for each chip on a night-by-night basis) 
and fixed. The WFC is known to have non-linearities in all 
its chips, and these were corrected for using the following 
equations (Mike Irwin, private communication): 



CC Dl,true 
CC 'i?2,true 
CC D^ t true 
CC Dijtrue 



= x - 2.5 x Kn 6 x 2 + 1.2 x 10~ n :r 3 



(2) 



1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

2 The authors acknowledge the use of the following software pro- 
vided by the UK Starlink Project: CCDPACK, CONVERT & 
KAPPA. Starlink is run by CCLRC on behalf of PPARC. 



Note that CCDi, t rue are the corrected counts (ADU) 
and that x are the raw, (bias-subtracted) values. 

Twilight flats were used as an initial correction for the 
images. However, sky fiats were needed for the full correc- 
tion and these were constructed from our object images. As 
can be seen in figure each final mosaic is made up of two 
pointings, one with the BCG centred on chip 4 (this is de- 
fined as pointing 1) of the WFC, the other centred on chip 
2 (pointing 2). Therefore, all chip 4 frames from pointing 1 
were automatically rejected from the construction process, 
as were all chip 2 frames from pointing 2. Every remaining 
frame was then inspected and rejected if it contained large 
bright sources or other defects. In order to keep signal-to- 
noise the same across the mosaic, the same number of frames 
was used to make the master sky flat for each chip. 

The master sky flats were constructed using a procedure 
based on that of Morrison et al. (1997) and Feldmeier et al. 
(2002). First, SExtractor (Bertin & Arnouts 1996) was used 
to pick out sources down to a threshold of la above the back- 
ground and create an object mask. The IRAF task mimstat 
was next used with this mask to calculate the modal value 
for each accepted frame. The frames were then normalised 
by their modes and median combined with IRAF's imcom- 
bine, using a sigma-clipping (±2<r) algorithm. Each of the 
individual sky frames was finally corrected with this mas- 
ter flat, and their modes were recalculated. The cycle was 
then restarted, with the original accepted frames now being 
normalised by their recalculated modes. This cycle of correc- 
tions successively reduces the distribution of sky values in 
the image and allows a more accurate determination of the 
modal value of the image. 4 such cycles were sufficient as by 
that point the percentage difference between modal values 
was of the order of 10 -8 . All object frames were divided by 
these final flats. 

Before combining the frames to make the final two- 
pointing mosaic, a correction was needed to account for the 
nonlinear 'pincushion' distortion inherent to the WFC op- 
tical system. Calibrated co-ordinate information specific to 
the WFC is available in the Starlink ASTIMP task in the 
CCDPACK suite of software. CCDPACK tasks were also 
used for the registering and mosaicing of the images, and 
for the level-matching between chips which was necessary 
due both to the different gain of each WFC chip, and for 
the background differences between the two pointings. 

3.2 Astrometry 

For each mosaic, an initial World Coordinate System was 
fitted using ~ 10 stars identified from Digitized Sky Survey 
(DSS) plates through the ALADIN Java applet (Bonnarel 
et al. 2000). Tasks from the WCSTools package (Mink 2002) 
were then used to locate a much larger number of stars (typ- 
ically ~ 100) from the USNO-A2.0 catalogue (Monet et al. 
1998). These positions were refined interactively and the fi- 
nal WCS was set by IRAF tasks using a fifth-order Legendre 
polynomial. Fitting the WCS allowed us to check that the 
BCGs being analysed in our images were the ones nearest 
the cluster X-ray peaks, an important point which has al- 
ready been mentioned above. The r.m.s. errors in the fitting, 
combined with the positional errors in the USNO-A2.0 cat- 
alogue (Deutsch 1998), give total uncertainities of the order 
of < 0".8 for image positions. 
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Figure 4. Masked images of the BCGs in A2256 (left) and Z4905 
(right). Both are r band images with the levels scaled logarithmi- 
cally. These images represent the extremes of the masking proce- 
dure, with Z4905 having 13.6% of the profile area removed from 
the isophote-fitting process and A2256 having 45.9% removed. 



Figure 3. Numbers of objects detected as a function of minimum 
area for the B-band image of A2029. The open circles show the 
distribution of detections in the positive image and the filled cir- 
cles show the distribution of detections in the negative image. The 
solid line is the difference between these and should be the distri- 
bution of real objects. The SExtractor minimum area parameter 
was chosen as the point where noise and object detections are 
equal. 



3.3 Masking 

Accurately modelling the BCG profile at very low surface 
brightnesses requires that contaminating sources are ade- 
quately masked. SExtractor was used to pick out objects for 
masking and there are two main parameters that need to 
be set for object detection: the detection threshold and the 
minimum number of contiguous pixels that constitute an ob- 
ject. We chose to mask down to a threshold of 0.8a above the 
background. However, at this level there is a danger of over- 
masking as it is easy to mistake noise spikes for real objects 
if the minimum object area is set too low. To determine an 
appropriate minimum area, for each image SExtractor was 
first run with a 0.8a threshold for a range of object areas 
from 5 to 20 pixels and the number of objects detected at 
each area setting was found. SExtractor was then run again 
in the same way on the negative version of the image. All 
objects detected in the negative image should be noise peaks 
and not real objects. If we assume that noise is distributed 
equally about the sky level, then the distribution of detec- 
tions in the positive image minus the distribution of detec- 
tions in the negative image will give the distribution of real 
objects as a function of minimum area (figure|5j. As a com- 
promise between masking too many noise spikes and not 
masking enough objects, the minimum detection area was 
finally set at the area at which the number of real objects 
was equal to the number of noise spikes being picked up. 
These detection limits correspond to a magnitude threshold 
(for point-like objects) of tub ~ 26.7 and m r i ~ 25.9. An 
examination of the number counts of objects in the masking 
catalogues shows that we are complete down to tub ~ 25.8 
and 77V ~ 25.0. 

SExtractor was then run to produce both an object cat- 
alogue and a segmentation image. Shape information from 
the output catalogue was used to produce an object image, 
with the major and minor axes of each object being scaled 



by its magnitude. The object and segmentation image were 
then combined to produce the final mask. The same mask 
was used for both B and r' versions of each image to ensure 
consistent colour measurements, with the mask being based 
on the r' image as it is deeper than the B image. Similarly 
to Krick et al. (2006), for each cluster two additional masks 
were produced with the magnitude-scaled major and minor 
axes shrunk by 10% and grown by 10%, in order to investi- 
gate the effect mask size has on the extracted BCG profile. 
Table 3 gives the percentage of each image masked. Typi- 
cally the masking procedure removes ~ 25% of the pixels 
from the images around each BCG, as listed in Table 2. The 
third column in the table gives the proportion of the area 
masked around the BCG out to the maximum radius of the 
software-generated BCG model. Figure 2] shows the 'best' 
and 'worst' cases from this column: The r' band image of 
the BCG in A2256 has ~ 46% of its area masked whereas 
only ~ 14% of the area around the BCG in the r' band 
image of Z4905 is masked. 

3.4 Sky subtraction 

Accurate sky subtraction is crucial when attempting to mea- 
sure surface brightnesses down to the level we require. We 
used an iterative procedure to remove the sky, which acted 
to improve our sky model with each cycle: 

(1) A large region around and including the BCG was 
masked out to a radius of 0.5Mpc from the BCG centre. 
This value was chosen as we did not expect to find any 
ICL beyond this at the depth of our images (Gonzalez et 
al. 2005, Zibetti et al. 2005). This was combined with the 
mask produced through the technique described above and 
used to reject points from a surface-fitting task (Starlink's 
SURFIT). The extent of the background sampled by our 
images (~ 41' x 41') deems it necessary to allow for large- 
scale gradients in the sky level and so a planar surface was 
fit to the background. 

(2) The IRAF ellipse task (Busko 1996, Jedrzejweski 
1987) was used to model and remove any large bright stars 
and galaxies (including the BCG) from the image. The el- 
lipse ellipticity parameter is normally required to be between 
0.05 and 1 but the parameter set was edited to allow it to 
be as low as 0.001 (the ellipse-fitting algorithm diverges at 
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Table 3. Extent of masking. Column 2 gives the percentage of 
the total area of the image that has been masked and the values 
in column 3 give the percentage of the profile area that has been 
masked in the B and r' bands, respectively. These values are 
different as profiles often extend further in r' . 
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e = 0) and fixed at this value to accurately model the cir- 
cular stellar isophotes. 

(3) The surface previously subtracted was then added 
back to the source-removed image and a new planar sky was 
calculated, as a better estimate of the actual background is 
made possible by the removal of bright sources. 



This cycle was run 3 times for each image. For the first 
of these cycles, an additional step was included between (2) 
and (3) in which a new mask was produced. This mask cov- 
ered the diffraction spikes and saturated regions of the bright 
stars removed in step (2) and allowed better stellar models 
in subsequent cycles. 

In order to check our large-scale flat-fielding and sky 
subtraction errors, we followed the approach outlined in 
Feldmeier et al. (2002). We binned up our final reduced, 
sky-subtracted images into 50 x 50 pixel sections and calcu- 
lated the median value of each bin. Each bin had to have at 
least half its pixels unmasked for its value to be accepted, 
and histograms of these medians were then created. Two ex- 
amples are shown in figure[S] The widths of these histograms 
around zero are in effect a measure of the above errors, as the 
dispersion is due to flat-fielding errors and excess contam- 
inating light which we have failed to mask. The mean sky 
values for the images shown, which are of the clusters A1672 
and A2029, were 437.2 ADU and 1145.5 ADU respectively, 
which we estimate give us uncertainties of 0.16% (a = 0.69 
ADU) and 0.09% (a = 1.04 ADU) for residual flat-fielding 
and sky-subtraction errors. The uncertainties for the full 
sample are provided in table 2. 

3.5 Comparison with Star Subtraction 

Instead of masking all objects in an image, the usual ap- 
proach to taking care of contaminating sources in an image 
is to first subtract off the stars, and then mask out the re- 
maining unwanted features. We opted for full masking as it 
was a technically less complex approach and we judged it 
to be less likely to introduce errors in sky subtraction. An 
incorrect PSF model will have a detrimental effect on the 
sky model as stars may be under- or oversubtracted, but 
this is avoided by simply masking out the stars and reject- 
ing them from the fitting. The main driving factor, however, 
was the speed of the masking process compared to one in- 
volving star subtraction, which is an important considera- 
tion for our project due to both the large size and number 
of our images. 

To ensure that we were not sacrificing accuracy for 
speed, for a few test images we implemented a star subtrac- 
tion algorithm. The basis for this was to produce a good PSF 
model that would account not only for the central regions 
of stars but also for the extended outer regions, as faint 
residual light is our main problem. As in Feldmeier et al. 
(2002) and Gonzalez, Zabludoff & Zaritsky (2005), we con- 
structed a large-radius PSF using the cores of unsaturated 
stars and the wings of brighter, saturated stars as these con- 
tain high signal-to-noise. This model was used to create a 
star-subtracted image which was used as the initial input 
for the iterative masking process described in §3.3. Compar- 
isons of the final ID profiles calculated by ellipse using the 
star-subtraction and masking-only methods are shown in fig- 
urc|SJ for two of our tests, A2249 and A2029. It is clear that 
the profiles agree extremely well, even down to the lowest 
surface brightnesses we can reach. 

4 RESULTS 

For each of the observed BCGs, figure ?? shows a 300" x 300" 
image in r' . 
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Figure 5. Histograms of background values for the images of A1672 (B) (left) and A2029 (r') (right), calculated after flat- 
fielding and sky subtraction. The bins have widths of 0.25ADU and the overlaid Gaussian fits have an r.m.s. of 0.69ADU 
and 1.04ADU, respectively (corresponding to surface brightnesses of p,B = 29.34 and fi r = 28.51). These histograms allow 
us to quantify our large-scale errors, which need to be as low as possible to reliably detect faint ICL. The middle plots 
give the median value of the bins as a function of distance from the centre of each image, whilst the lower plots show the 
running median of lOOpixel radius intervals. They show that there appears to be no significant correlation of residual value 
with position. 



Radial profiles were calculated from the masked image 
in each pass-band for all target galaxies, using the mask 
produced from the r' band image. We used the masking- 
only method of removing foreground objects. As discussed 
in §3.4, this method is as reliable as subtracting stars using 
a PSF fit, and is also faster and more robust. As mentioned 



in §3.3, we used the IRAF ellipse task to make the measure- 
ments. The radii of the isophotes were chosen to follow a ge- 
ometric progression and the centroids were kept fixed. The 
resulting surface brightness profiles are plotted in figure ??, 
along with colour, ellipticity, and position angle profiles. To 
ensure consistency in colour measurements, colour profiles 
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Figure 6. Plots of the difference between the surface bright- 
ness profiles calculated using the star-subtraction and masking- 
only methods, for the BCGs in A2029 (top) and A2249 (bot- 
tom). See main text for an explanation of the differences in 
these approaches. As the profiles show, the simpler and less time- 
consuming masking-only approach produces profiles which agree 
very well with the star-subtraction method, even at the faint ends 
of the profiles, with no appreciable systematic difference in mag- 
nitude with radius. The error bars are based on the uncertainty 
in isophotal intensity. 



were calculated using the geometry information for the r' 
band isophotes for both the r' and B images. 

Error limits for the surface brightness and colour pro- 
files come from combining large-scale errors (§3.4), an es- 
timate of star/galaxy light missed from the masking proce- 
dure (using the 'grown' and 'shrunk' masks) , and the scatter 
of intensities along the isophotes. In virtually all cases, the 
last two sources of uncertainty are negligible, even at the 
faint ends of the profiles. Instead, the uncertainty is domi- 
nated by the effects of flat-field inhomogeneities and imper- 
fect sky subtraction. However, the errors shown in the pro- 
files and table 2 will be an overestimate of the uncertainty 
as we are considering variations and missed light across the 
entire image not just the area covered by the host cluster of 
the BCG. 
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Figure 7. Offset between the BCG centre and the X-ray emission 
peak of its host cluster (from the BCS). The inner dashed circle 
has a radius of 50kpc and the outer is lOOkpc. The size of each 
marker is the error in position associated with the calculated off- 
set. The mean offset between the peaks and centres is 129.0kpc, 
with an r.m.s. of f22.0kpc. 



5 DISCUSSION 

We have presented basic data in the form of surface- 
brightess, ellipticity and position angle profiles for a large 
sample of BCGs. In a further series of papers we intend to 
perform a variety of analyses on our data set, including an 
examination of the morphological properties of the BCGs, 
the investigation of substructure, and surface brightness pro- 
file fitting. 

For our 49 clusters we find that in 22 cases the cen- 
tre of the BCG is less than lOOkpc from the centre of the 
X-ray isophotes taken from the BCS survey (Ebeling et al. 
1998, 2000). Cypriano et al. (2004) found a similar align- 
ment for 17 of the 22 clusters they examined for a weak 
lensing study. We note that the cluster X-ray positions ob- 
tained from the NORAS catalogue (Bohringer et al. 2000) 
have an RMS scatter of 95kpc compared to the BCS po- 
sitions, so this scatter is largely related to the difficulty of 
defining the isophotal centre. In two cases there is a dis- 
crepancy in excess of 0.5Mpc between the purported X-ray 
centre and the location of the BCG. Closer examination of 
these cases reveals that both these objects have significant 
levels of substructure with obvious large sub-clumps. Such a 
close agreement between the centre of the mass distribution 
and the location of the BCG indicates a cuspy mass profile 
with a small or negligible core, in agreement with modern 
simulations of large clusters (Navarro, Frenk & White 1996; 
Moore et al. 1998; Power et al. 2003) which readily produce 
such profiles. It also supports the contention of Allen (1998) 
that the discrepancy between X-ray and gravitational lens- 
ing mass estimates for clusters of galaxies was due to an 
overestimation of the core size in so-called non-cooling flow 
clusters that was not present in the underlying mass dis- 
tribution. Such a conclusion had already been reached by 
Smail et al. (1995) using strong lensing measurements. 

It is thought that BCGs form early in the history of 
their host clusters (Merritt 1985, Tremaine 1990) and that 
in general they will have grown little since then. Work by 
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Figure 8. Surface brightness and Petrosian R p (= 1 profiles for the BCGs in A2124 and A2249. Plotting R p to examine the behaviour 
of the BCG light profile means that determining the limits of the Sersic fit becomes much less subjective, which is particularly important 
when the BCG appears to have a halo, as is the case for A2249. Note that the innermost dashed lines denote the seeing radii and the 
outer dashed line in the plot of A2249 is the radius outside which data is excluded from the Sersic fit. 



Conselice, Bershady and Jangren (2000) suggests that quan- 
tification of observed asymmetries coupled with colour in- 
formation may allow us to determine whether a galaxy is 
involved in an interaction and therefore allow us to see how 
often events that contribute to the growth of the BCG ac- 
tually occur at low redshifts. 

Substructure within the cD envelope/ICL will also en- 
able us to study the history of interactions. Large tidal de- 
bris arcs have been found in the Coma cluster (Trentham 
& Mobasher 1998; Gregg & West 1998) and the Centau- 
rus cluster (Calcaneo-Roldan et al. 2000). Feldmeier et al. 
(2002) also report a tidal debris plume in MKW7 and fur- 
ther features in A1914 (Feldmeier et al. 2004). It has been 
proposed that these features arise from tidal interactions be- 
tween cluster galaxies and the cluster potential (Moore et al. 
1996) and such arcs have been seen in simulations following 
the evolution of intracluster stars (Willman et al. 2004). 

It has long been known that BCGs have shallower sur- 
face brightness profiles than the de Vaucouleurs r 1 ' 4 profile 
(de Vaucouleurs 1948), and it is this excess light over the r 1 ' 4 
profile which is usually termed the cD envelope. Beginning 
particularly with the work of Caon, Capaccioli & D'Onofrio 
(1993), many authors (see Graham & Driver, 2005, for an 
extensive list) are now favouring the r 1//n , or Sersic (Sersic 
1968), profile as the universality of the r 1 ^ 4 is questionable. 
In addition, the exponent n (often known as the shape pa- 
rameter) , rather than being simply an extra parameter with 
no physical basis, has been found to correlate with various 
other observable properties in a way that is not explained 
by parameter coupling, thus offering a much deeper insight 
into bulge-type systems than the r 1 ^ 4 law. The value of n 
correlates with effective radius r e and the total luminos- 
ity of the system, such that more luminous galaxies have a 
larger value of n (e.g. Caon, Capaccioli & D'Onofrio 1993); 
Graham, Trujillo & Caon (2001) show that log(n) corre- 
lates with the the galaxy's central velocity dispersion; and 
n, effective radius and central surface brightness are tightly 



distributed about a plane (the photometric plane) for ellip- 
ticals (Khosroshahi et al. 2000; Graham 2002; La Barbera 
et al. 2004; La Barbera et al. 2005; Ravikumar et al. 2006). 

However, Sersic fitting suffers from the following prob- 
lem: given the presence of an envelope, the range of radii over 
which the model should be used to fit the surface brightness 
profile is difficult to determine. Although this problem also 
occurs with the r 1 ' 4 law (Schombert 1986), it is made worse 
in the Sersic case because of the extra degree of freedom. An 
interesting solution to this that we have found makes use of 
the Petrosian index, r\ (Petrosian 1976). This is the ratio of 
the average intensity within some radius (R) to the intensity 
at that radius: 



V (R) = 



2f*I{R')R'dR' 



R 2 I{R) 



(3) 



Figure |H] shows surface brightness profiles of the BCGs 
in A2124 and A2249, and their corresponding plots of R p 
(— I/77) as a function of radius. The monotonic decrease 
in the R p profile of A2124 implies that there is very little 
excess light found here, i.e. this galaxy does not appear to 
have an extended halo, and this can be seen in its surface 
brightness profile which is well fit at all radii by a Sersic 
law. On the contrary, the other BCG shows an increase in 
-Rp at a radius of ~ 30". The point at which R p increases 
can be translated as the point where we begin to see an 
envelope, and where the Sersic fit to the surface brightness 
profile becomes unreliable. Note that the innermost dashed 
lines denote the seeing radii and the outer dashed line in the 
plot of A2249 is the radius outside which data is excluded 
from the Sersic fit. However, we note that although we do see 
a correspondence in these two examples, a rigorous analysis 
is required to put this method onto a reliable footing, and 
we intend to carry this out in future work. 

35 of the 49 BCGs show a very clear increase in el- 
lipticity with radius, and this may mark the change from 
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the stellar population of the BCG to that of the intracluster 
medium. Basilakos, Plionis & Maddox (2000) found that the 
distribution of ellipticities of clusters in the APM survey has 
a peak at e = 0.46. This agrees well with the values for the 
outermost isophotes of many of our BCGs, suggesting that 
perhaps the outer regions are tracing the cluster potential. 

As we have already mentioned, one of the major ad- 
vantages of our data set is that we have imaging in both 
B and r' bands, allowing colour gradients to be measured. 
Schombert (1988) considered two-component fits in order to 
derive properties for the cD envelope, and this was taken a 
step further by Gonzalez et al. (2005), who used full two- 
dimensional two-component fits. It will be extremely inter- 
esting to combine these approaches with colour analysis to 
understand the properties of the BCG and its envelope. 
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